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Abstract 

We use the polynomial formulation of the holomorphic anomaly equations governing per¬ 
turbative topological string theory to derive the free energies in a scaling limit to all orders in 
perturbation theory for any Calabi-Yau threefold. The partition function in this limit satisfies 
an Airy differential equation in a rescaled topological string coupling. One of the two solutions 
of this equation gives the perturbative expansion and the other solution provides geometric hints 
of the non-perturbative structure of topological string theory. Both solutions can be expanded 
naturally around strong coupling. 
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1 Introduction 


The perturbative expansion of string theories is asymptotic [1, 2] which raises questions about 
the non-perturbative completion and definition of these theories. It has been very fruitful to 
address these questions within topological string theory which can be connected to Chern- 
Simons theory and matrix models, see the excellent review [3] and references therein. The 
lessons obtained from this study are expected to give general insights into string perturbation 
theory. 

Topological string theory is based on the nonlinear sigma model which maps the string 
world-sheet into a target Calabi-Yau (CY) threefold. The topological string partition function 
is given by a perturbative expansion in the topological string coupling A, summing the free 
energies jr{g) 

over the world-sheet genera g: 

oo 

Ztop = exp E 
9=0 

The perturbative free energies satisfy the holomorphic anomaly equations of Bershadsky, 
Cecotti, Ooguri and Vafa (BCOV) [4, 5] which are recursive differential equations. These were 
solved using Feynman diagrams [5], which give the form: 

E^-m+ /'»>. 0.2) 

redMg 

where is the Deligne-Mumford compactification of the moduli space of Riemann surfaces 
of arithmetic genus g. The stratification of this moduli space can be captured by the decorated 
dual graphs T of the Riemann surfaces. The first summand in the above equation is a summation 
over the decorated graphs corresponding to the degenerate Riemann surfaces. The weight o;(r) 
is given by the Feynman rules: for example, an edge in the dual graph corresponds to a node in 
the degeneration and gives a propagator. The contribution to from the smooth Riemann 
surfaces of arithmetic genus g is the holomorphic ambiguity which is fixed by boundary 
conditions. The factors l/|Aut(F)| in the above sum are universal constants coming from the 
structure of M.g and are independent of the geometric data of the target CY threefold whose 
information is encoded in the weights a;(F) obtained from the Feynman rules. 

The Feynman diagrams become intractable at higher genns since their number grows rapidly. 
It was shown that this simplifies signihcantly due to a polynomial structure of the higher genus 
free energies in finitely many generators [6, 7]. The polynomial structure is non-trivial and 
stems from further decomposing the vertices of the Feynman diagrams at higher genus into 
simple expressions as well as from a differential ring structure of the generators. This leads to 
the much milder polynomial growth of the number of terms at each genus. In this paper we 
use this polynomial structure to determine at all genera the pieces of the free energies coming 
from certain monomials at every genus g which include contributions from the most degenerate 
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Riemann surfaces in dA4g. To select these terms in the full partition function we rescale 
the polynomial generators as well as the topological string coupling such that only these terms 
survive a scaling limit. We derive a modified Bessel differential equation in a rescaled topological 
string coupling whose solution assembles these monomials at all genera. This solution allows us 
furthermore to expand this piece of the partition function around strong coupling A —)• oo. The 
second linearly independent solution of the differential equation is non-perturbative at A = 0 
and provides geometric hints of non-perturbative topological strings. 

Furthermore, a change of variables allows us to transform the differential equation into an 
Airy differential equation. The latter is the hallmark of two dimensional topological gravity 
studied using matrix models [ 8 , 9] and appears in many discussions of non-perturbative phe¬ 
nomena [3] and especially also in the recent definition of non-perturbative topological strings 
on some non-compact CY manifolds given in Ref. [10]. It is surprising that the results of the 
present work are universal for topological strings on any CY geometry which is subject to the 
holomorphic anomaly equations even if there is no manifest relation to Chern Simons theory or 
matrix models such as compact CY threefolds. 

The polynomial structure of the higher genus free energies was already used in Ref. [ 6 ] 
to determine in principle the coefficients of certain monomials in the generators at all genera 
for the quintic. In a similar context, where the polynomial generators are realized as quasi 
modular forms for a non-compact CY given by the canonical bundle over x F^, this has 
been addressed in Ref. [II]. The mirror curve of the latter geometry is related to the to the 
matrix model description of ABJM theory [12] and has been used in Ref. [13] to obtain strong 
coupling results. In Ref. [14] the all genus free energies for the matrix model are summed up in 
a certain limit in the moduli space and the coefficients are related to the Airy function which 
is very close in spirit to our work.^ 

The structure of this note is as follows. We review the polynomial structure of topological 
strings in Sec. 2. In Sec. 3 we derive a differential equation which determines the partition 
function in a limit to all genera and discuss its transformation into the Airy equation as well as 
the strong coupling expansion of its solutions. We conclude with discussions in Section 4. 

2 Polynomial structure of topological strings 

In this section we review the polynomial formulation [ 6 , 7] of the holomorphic anomaly equations 
[5]. Let Ai be the moduli space of a Calabi-Yau threefold, which can be the moduli space of 
complexified Kahler structures of a CY threefold Y or the moduli space of complex structures 
of its mirror X. In the following we will use local complex coordinates z*,i = l,...,n = dim A4 
and restrict to a one-dimensional slice of the moduli space by choosing a local coordinate z = z* 
such that C*** defined below is non-zero and we set z* = 0, z 7 ^ *. We use dz ■= -§^,82 ■= -§=■ 

^We would like to thank Marcos Marino for pointing out this reference to us and its relation to our results. 


3 



The function e~^ gives a Hermitian metric on a line bundle T —)• At with connection 
and provides the Kahler potential for the Weil-Petersson metric on Al, whose components and 
Levi-Civita connection are given by Gzz '■= dzdgK , = G^^dzGzz ■ The holomorphic Yukawa 

couplings are: Gzzz £ T (g) Sym^T*Al) , the curvature is expressed as : 

- i?,/, = [5„ DzYz = d,Kz = 2Gzz- - GzzzCY, (2.1) 

where Dz denotes the covariant derivative and ■= e^^G^^G^^Gzzz- This data defines a 
special Kahler manifold [15, 5]. 

We further introduce the objects which are non-holomorphic sections of C~^ (g) 

Sym^TAf with m = 2 , 1 , 0 , respectively, and give local potentials for the non-holomorphic 
Yukawa couplings: 


= c; 


d-zS^ = Gz-zS^ 


d-zS = Gz-zS^. 


( 2 . 2 ) 


These are the propagators of the Feynman rules derived for in Ref. [5]. 

The topological string amplitudes at genus g with n insertions are defined in Ref. [5] 
as non-holomorphic sections of the line bundles over Af. These are only non-vanishing 

for {2g — 2 + n) >0. They are related recursively in n by DzI-^YIi = ^n \ as well as in g by 
the holomorphic anomaly equation for 5 : = 1 [4] 

d-z^Y^ = \CzzzCt + (1 - y^)Gz-z , (2.3) 

where x is the Euler character of the CY threefold Y. As well as for g >2 [5]: 






DzF^^^DzF^^-^^ +DzDzF^3-^'> 


(2.4) 


It was shown in Refs. [ 6 , 7] that for any CY threefold Fn'^ is a polynomial of degree 2>g — 
3 -|- n in the generators , S^, S, Kz where degrees 1, 2,3,1 were assigned to these generators 
respectively. The purely holomorphic part of the construction as well as the coefficients of the 
monomials are rational functions in the algebraic moduli. The recursive proof of this relies on 
the differential ring structure of the generators. For the purpose of this work we will only need 
the following [7]: 

DzS^^ = 2S^ - GzzziSn^ + hf , (2.5) 

where denotes a holomorphic function which is fixed by a choice of satisfying Eq. (2.2). 
The expression for the curvature (2.1) can be integrated to: 


T 


2 

22 


2Kz - Gz 


qZZ I 2 


( 2 . 6 ) 


with a holomorphic function depending on the choice of We use this to write out the 
following: 


DzGzzz = dzGzzz - 3sl,Gzzz - ^KzGzzz + 2,GtzzS" 


(2.7) 
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The holomorphic anomaly equations split into two equations [7]: 
dT^9) 

~d^ 

0 = 

3 All genus differential equation 

We determine in the following the all genus coefficients of particular monomials appearing in the 
polynomial formulation of the free energies To this end we derive a differential equation 

in a rescaled topological string coupling for the partition function which can be transformed 
into an Airy equation. 

3.1 Scaling limit 

In the polynomial expression of we consider the highest degree term in the generator 
which is a monomial of the form: , where f{z) is a rational function of the modulus 

z. From the Feynman diagram rules of Ref. [5] and the polynomial structure reviewed in Sec. 2 
we know that it is of the form: f{z) = ag Ug G Q. 

The corresponding set of graphs in the Feynman diagrams includes in particular those special 
graphs F which correspond to the most degenerate Riemann surfaces of arithmetic genus g and 
of geometric genus 0. We denote this set by Adg, cubic- A Riemann surface in this set is obtained 
by gluing 2g — 2 genus zero Riemann surfaces, each one has three markings, along the markings 
pairwise. The dual graph is obtained by gluing the cubic vertices along the half-edges. The 
arithmetic genus is then the number of loops or equivalently the first Betti number of the dual 
graph. Among the configurations in Mg, these dual graphs have the largest possible number of 
loops. The monomials we focus on furthermore receive contributions from further decomposing 
the vertices which are given by higher genus amplitudes with insertions. 

We introduce the total free energy, omitting g = 0,1, 

OO 

5", 5, Kz; z,X) = Y^ 5, AT,; z ). (3.1) 

9=2 

To select the terms Og Czi ^'^from W we rescale the generators as well as the 
topological string coupling A with e in the following way: , A = j, and define: 

OO 

T,{\s) = \YmT{S^\S\S,Kz,z,\) = Y,ag\l^~'" , (3.2) 

9=2 

where we defined the rescaled coupling A^ = A^ . 


9-1 


^ h=l 

dF^9) 

+ -5 + -5 


2 


dKz 


dS 


dS^ 


( 2 . 8 ) 
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3.2 Differential equation for all genus free energy 

We will now use the holomorphic anomaly equations in their polynomial form (2.8) to derive 
an equation governing the coefficients {og} defined in Eq. (3.2). 

Proposition 3.1. The all genus free energy in the scaling limit iFg satisfies the differential 
equation: 

elFs + + 2 (l - + ^ = 0 , . (3.3) 

Proof. We first introduce the notation: 

+ l-o.t., g>2, (3.4) 

where we denote by l.o.t. all monomials which vanish in the partition function in the limit 
described above. We obtain further from Eqs. (2.5, 2.7): 

= ( 3 ^ _ ^ _ ^3 5 ) 

At <7 = 1, the integration of Eq. (2.3) becomes: 


:= + l.o.t.. 


(3.6) 


At g = 2, multiplying both sides of Eq. (2.8) by fixes 02 = ^. For g > 3 the L.H.S. becomes: 


,dF^9) 

~d^ 


= a5(39-3)C2.-2(S“)='»-» + l.o.t. 


(3.7) 


and the R.H.S. is: 


^ (^^ahag-h{3h - 3)(3ff - 3h - 3) + ag.ffSg - 6){3g - 3)^ + l.o.t.. (3.8) 


Eq. (3.3) is obtained from the summation: 


^ ffiTrig) QZZ 


^g-l 


9=2 




- 1 ) 


(3.9) 


9=2 


\h=l 


□ 


3.3 Modified Bessel equation 

The equation for the partition function Ztop,s = exp/), in the scaling limit becomes: 
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Proposition 3.2. 2top,s{^s) satisfies the following differential equation: 

3Ag 3A^ Q 



1 

3A^ 


(3.10) 


This is the modified Bessel differential equation in terms of the variable ^ and the general 
solution in terms of the modified Bessel functions is given by: 


Z’top^S - 


A. 


"■h ( slf ) +“"■''■'11 ~ 


1 (sll)) ■ 


(3.11) 


Now we discuss the asymptotic behavior of the modified Bessel functions (see e.g., Ref. [16]). 
If Cl = 0 the series around A^ = 0 coming from ilii/s is trivial. If ci ff 0, the asymptotic series 
expansion around As = 0 is independent of C 2 . This gives the sequence {ag}g >2 up to a constant: 


J^s-{ 3^2 In As 


(1 .X 5 ,2 5 ,4 1105,6 

+ ( ^ + In As + -A, + -A, + —A, + 


(3.12) 


The first two terms come from the prefactor in Eq. (3.11) and can be regarded as the contribution 
of genus zero and genus one free energies to J^s while the rest from the modified Bessel function 
part. We fix ci = 1 and parameterize the general solution with an additional parameter C, by: 


^top,S - 


As 


h 

3 


cec. 


(3.13) 


In particular f does not affect the perturbative expansion around As = 0 but becomes relevant 
away from this locus. It can be interpreted as giving a non-perturbative correction to the 
perturbative series of topological string partition function. 

3.4 Strong coupling expansion 

Eq. (3.10) has two apparent singularities: the one at As = 0 is an irregular singularity, while 
the one at As = oo is a regular singularity. Expanding Ztop,s near As = oo, we get 


■'top,s 


^ ^/3^) ( 3 A 2 ) + ( 3 A 2 


TT 




00 z' 1 1 X 2 n 

- ^ 12 SA? 1 




1 27r ^ — K 

+ 63—Ce 

v3 


n =0 + 3) 

00 /I 1 \ 2 n 
(23A^) 


n=o’"!r(n+3) 


Up to irrelevant factors, one has near A = 00 the following: 

2 _ TT ^ ^ 

- - In As + 6 3 —As ^ + C)(A^ ^). 


(3.14) 


(3.15) 
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3.5 Airy equation 


2 1 

We can further make the change of variables: z = ( 2 A^)“ 3 ,t) = 2~3e^XiZtop,s- Eq- (3.10) 
then becomes the Airy equation; 

[d‘^ — z) v(z) = 0 . (3.16) 

This offers a more geometric picture of the non-analytic behavior of the modified Bessel func¬ 
tions. Applying the Laplace transform, one obtains: 


v{z) = — [ > 
^ ’ 27Ti 




(3.17) 


The integral contour 7 on the tc-plane has to be chosen such that the integrand vanishes at the 
boundary. There are essentially three choices for them given by 7 j,i = 1,2,3, satisfying the 
homology relation 71 -|- 72 + 73 ~ 0 as depicted in Fig. 1. 



Figure 1: The contours 71,72 and 73 on the tc—plane. 


The Airy function Ai( 2 ;) as a solution to this equation corresponds to the modified Bessel 
function Ki and is given by: 


Ai(z) = 


7 r \/3 3 


2 3 

3^" 


1 

27ri 


i j—-f2 






(3.18) 


It has exponential decay as z —)■ 00 , arg z = 0 and is oscillating as z —)■ 00 , arg z = tt. The other 
independent solution is usually taken to be the function 


Bi(z) 




+ 7_i 

3 



1 

27rz 


,^zwe-\w^dw. 


'172+2173 


(3.19) 


It has exponential growth as z —>■ 00 , arg z = 0 and is oscillating as z —>■ 00 , arg z = n. The same 
integral would have different asymptotic series expansions as the phase of z changes, exhibiting 
Stokes phenomena as discussed in e.g.. Refs. [3, 17]. 







The partition function can be written as: 

Ztop,s = 2h~^sx;K • (3-20) 

The integral over —72 gives Ai( 2 ;) and thus the purely non-perturbative part of Ztop,s, while the 
integrals over the other two cycles 71 , 73 are equivalent modulo the non-perturbative expressions 
and contribute to the perturbative part of Ztop,s- Analytic continuation and exploring the non- 
perturbative content can thus be realized via moving the integral contour as discussed in many 
other contexts, such as in Ref. [17]. 

4 Conclusion and discussion 

In this work we determined certain terms of the topological string free energies to all orders in 
perturbation theory universally for any CY threefold. Using the polynomial formulation of the 
holomorphic anomaly equations we derived an Airy equation in a rescaled topological string 
coupling which has a solution encoding these terms. This solution admits a strong coupling 
expansion which should be the analog of an expansion considered in Ref. [13] for ABJM theory 
[12]. It would be exciting to develop a similar interpretation which can be attached to any CY 
background. 

A second linearly independent solution of the Airy equation does not contribute to the 
perturbative expansion around = 0 but gives contributions away from this value. This is a 
manifestation of a non-perturbative ambiguity attached to differential equations as it appears 
in many other contexts, see Ref. [3]. This may offer geometric hints of the non-perturbative 
structure of topological string theory. Indeed the Airy function appears in the recent work [10] 
defining non-perturbative topological strings for some non-compact CY threefolds. 

Differential equations in the parameters of a theory often suggest an underlying geometric 
interpretation and perhaps a variation problem associated to it. Our results may provide the first 
steps in this direction. It is perhaps suggestive to think about the topological string coupling A 
as a complex variable, paralleling the discussion of analytic continuation of Chern-Simons theory 
[17]. Furthermore, the non-perturbative structure of topological strings was also addressed using 
methods of trans-series and resurgence of differential equations, see e.g.. Ref. [18] and references 
therein. These methods were applied to the master anomaly equation which is not a differential 
equation in the topological string coupling. Our results show that a differential equation in the 
coupling can be deduced, giving further motivation to pursue this line of research. 

The limit described in the paper seems to select only the cubic graphs or equivalently the 
geometric genus zero contribution to the topological string partition function. The similarities 
to Refs. [ 8 , 9, 19] suggest to further study an associated integrable hierarchy structure of the 
partition function. It is natural to expect, from examining the Feynman diagrams, that the 
lower degree terms in the polynomial generator should correspond to Riemann surfaces with 
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higher geometric genera. In this sense Ztop,s is really the classical part with respect to the loop 
expansion in the parameter e assigned to the geometric genus, and the full partition function 
Ztop is its quantization. This parameter plays an analogous role as the equivariant parameter 
in localization. The interpretation of this parameter in the A-model would have the effect 
of reducing the equivariant localization to the stratum corresponding to the most degenerate 
conhgurations in the moduli space of curves. It would be furthermore interesting to understand 
a more direct physical meaning of the parameter e. These questions will be addressed elsewhere. 
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